Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1140 1823 (N56E17 L= 36)  1644    242        1 (1044/1772-1228/1853)  56.3 E24-E08 N48-N63
      2 1787 1127 (N21W52 L=104)  1132    157        3 (1739/1073-1830/1196) -54.5 W46-W57 N17-N24
      3 1701 1154 (N20W43 L= 95)   386     74        1 (1665/1126-1729/1181) -18.4 W39-W45 N18-N21
      4  222 1134 (S10E60 L=352)   766    139        2 ( 181/1072- 270/1178)  45.0 E64-E54 S14-S07
      5 1650 1113 (N17W39 L= 91)   866    198        3 (1589/1024-1698/1173) -28.6 W33-W44 N12-N19
      6 1565 1098 (N14W32 L= 85)   277     67        2 (1547/1069-1579/1134)  78.7 W30-W33 N11-N16
      7 1746 1045 (N15W48 L=101)   624    115        3 (1718/ 986-1783/1085)  66.8 W45-W52 N10-N17
      8 1215  634 (S21W21 L= 73)   272    115        0 (1170/ 602-1263/ 669)  56.3 W18-W22 S23-S17
      9 1699  512 (S17W64 L=117)   651     80        0 (1674/ 478-1735/ 535)  53.1 W61-W68 S19-S15
     10 1466  446 (S27W46 L= 99)   908    164        3 (1398/ 426-1538/ 482)  41.2 W41-W50 S29-S22
     11 1280  385 (S35W36 L= 88)   374     90        1 (1239/ 369-1324/ 405)  38.7 W33-W38 S36-S32
     12 1151  240 (S50W38 L= 91)   479     51        0 (1118/ 215-1177/ 260)   0.0 W33-W42 S52-S48
     13 1311  196 (S49W61 L=114)   789     92        2 (1273/ 164-1338/ 218)  -5.2 W52-W74 S50-S46
  

    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