Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1411 1562 (N32W39 L= 32)   561    183        2 (1330/1541-1476/1592)  -9.5 W32-W44 N30-N33
      2 1109 1441 (N30W13 L=  7)   848    179        1 (1044/1387-1208/1479)   8.7 W07-W20 N27-N31
      3 1672 1356 (N14W51 L= 45)  1348    183        2 (1599/1292-1748/1396) -45.0 W45-W55 N08-N17
      4 1146  956 (N00W06 L=360)  1739    459        7 (1047/ 774-1198/1171) -90.0 W03-W08 S10-N14
      5  526 1076 (N16E31 L=323)   504    152        0 ( 465/1047- 593/1114) -29.1 E35-E26 N13-N19
      6 1426  918 (S07W23 L= 17)   535    113        2 (1387/ 886-1475/ 954) -39.8 W20-W26 S08-S03
      7 1319  798 (S12W14 L=  8)   386     94        0 (1298/ 745-1342/ 851)  66.8 W12-W16 S15-S08
      8 1327  589 (S26W13 L=  7)   623    133        1 (1299/ 523-1359/ 649)  54.5 W10-W15 S29-S22
      9  973  345 (S37E17 L=337)  1200    291        1 ( 891/ 280-1024/ 453)  40.6 E25-E11 S42-S29
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y